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Abstract 

In this work, the transition matrix elements for inelastic electron-electron scattering are inves- 
tigated. The angular part is given by spherical harmonics. For the weighted radial wave function 
overlap, analytic expressions are derived in the Slater-type and the hydrogen-like orbital models. 
These expressions are shown to be composed of a finite sum of polynomials and elementary trigono- 
metric functions. Hence, they are easy to use, require little computation time, and are significantly 
more accurate than commonly used approximations. 
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I. INTRODUCTION 



Ever since Geiger and Marsden shot a particles on a gold foiP and Rutherford subse 



quently used their results to confirm that atoms consist of a very small nucleus surrounded 
by an electron cloucP, scientists have used scattering effects to determine the properties of 
otherwise invisible or inaccessible objects. 

In particular, quantum systems are usually investigated by means of scattering experi- 
ments. The central quantity in these quantum mechanical scattering systems is the matrix 
element of the transition operator V , 

{f\{F\V\I)\^. (1) 

which describes the amplitude for the transition from a product state |/) \i) to another 
product state \F) \ f) in first order Born approximation. Here, we distinguish between the 
states of the probe (lower-case letters) and the target (upper-case letters).'^ 

In this work, we limit ourselves to the treatment of the electronic transition of a single, 
isolated atom, and the effect it has on the probe triggering that transition. Assuming a 
Coulomb-like transition potential, 

(r'l (R'lVlR) |r) oc , ^ , ^(r^ - y)5CR! - R) 
|r — R| 

(k'l (R'l^lR) |k) oc - R) (2) 

piQR 

with Q := k — k' and inserting appropriate identity operators, the transition matrix element 
can be written as 

(/I {F\V\I) \i) 

(/|k') (F|R') (k'l (R'|V|R) |k) (R|/) (k|.) 

(/|k') (F|R)e^Q-^(R|/) (kK). 



R,R',k.k' 
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Here, we assume that the probe electron was prepared as a plane wave, i.e., \i) = |k), and 
that the detector also measures plane waves, i.e., I/) = I k') .HI In that case, the transition 
matrix element assumes its commonly used form, 

i^^(F|R)e^Q-^(R|/), (4) 



or, by virtue of the partial wave or Rayleigh expansion, 

% E i'n"(Q/Q)* (m) Y,^{R/R)UQR) (R|/) , (5) 



where the Yj^ denote the spherical harmonics and the jx are the spherical Bessel functions 
of first kind. 

Up to this point, the treatment is exact. In order to evaluate the matrix element further, 
one needs to explicitly specify the functions (R|/) and (R|-F), however. For the electronic 
transitions treated here, two common approximations are Slater-type orbitals (STO)0S1 and 
hydrogen-like orbitals (HLO)^. In both models, the initial and final states are modeled as a 
spherical harmonic angular dependence and a (analytic) radial wavefunction ip: 

(R|/) = ^i{R)YriB./R) (R|F) = ^p(i?)F-'(R/i?). (6) 



With this, eq. |5] becomes 



oo 



(7) 



xIjf{R)3x{QR)^AR)R^<^R, 

where the angular integral evaluates to Wigner 3j symbol^^. Hence, we will only deal with 
the radial integral 

/■oo 

(jA(g)) := / ^F{R)jx{QR)iJi{R)R^(^R (8) 

Jo 

in this work, for which we will derive an analytical form in the STO and HLO models. 

The rationale behind this choice of orbitals is the fact that the initial state is usually a 
tightly bound (and hence strongly localized) core state that is described well by an atomic 
model. Since the matrix element can only be non-zero if both the initial and final states are 
non-zero, this effectively selects a portion of the final state that is close to the nucleus as 
well. 

Crystal and many-body effects, on the other hand, are often relatively small perturbations- 
especially in the commonly experimentally accessible regions. In addition, for elucidating 
the underlying fundamentals, an analytical treatment in an isolated atom model is usually 
advantageous. 

In many articles, eq. [s] is simplified further by taking the small angle limit (also known 
as dipole-approximation). 



in which case the integral over R in eq. [s] becomes only a weighting factor (the factor 
can be moved out of the integral). 

For increasing momentum transfer, a leading term Q'^ would increase boundlessly, giving 
rise to a similar effect as the ultraviolet catastrophe. To avoid that, an artificial cutoff Qc is 
sometimes introduced, which is only fighting the symptoms instead of the cause, however, 
and is not very elegant. In addition, is has recently been shown experimentally using electron 
energy loss spectrometry (EELS) that a dependence is an oversimplification even for 
Q < Qc and can lead to errors of the order of 25 9c!^m^. 



where n is an effective quantum number (which is not necessarily an integer). C, := {Z — s)/n 
is an effective nuclear charge, with the physical nuclear charge Z and a screening factor s. 
The normalization constant is given by 



II. SLATER-TYPE ORBITALS 



The radial part of STOs is given bj^^ 



(10) 




(11) 



Then, eq. [8] becomes 




(12) 



where we defined n := ni + np and C, = (C/ + Cf)/o^- As derived in appendix [A| this integral 
can be evaluated analytically, yielding 




(13) 
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TABLE I. Weighted radial wave function overlap in the STO model for monopole (A = 0), dipole 
(A = 1), and quadrupole (A = 2) transitions. For the definitions of the constants Nj,Np,nX refer 
to the text. 



_\ = Q ^'^^'^["'l, sin I n arctan 



A = 1 



NjNFr{n-l) 
Q2(C2+Q2)5 



C sin ( n arctan ( ^ ) ) — Qn cos ( n arctan ( § 



A = 2 QHe+Q^)"^ 



C{Q^ - 3nQ^ - n^Q^ + 3^) sin n arctan 



+ 3nC - n^Q'^ + 3C^) cos ( n arctan 



Q 



with the constants A^^,B^^ as defined in appendix [a| This is the first main result. 

In table [T| eq. [13] is evaluated for the most important transitions with A = 0, 1,2, corre- 
sponding to monopole, dipole, and quadrupole transitions, respectively. 

It should be noted that for Q — > oo, with 



sm 



cos 



a arctan 



a arctan 



sm 



cos 







/air 






It 


a7[\ 




fan 









(14) 



it is straight forward to derive that eq. [13] behaves as 
NiNf 



Q 



n+l 



A+1 



sin((n-A)|)^(-l)^(At,„ + i?y 



fc=0 



A + 1 



c 

Q 



cos 



((^ - ) ((2fc + ^ - ^)4n + {2k + n-X- l)Bl^) 



k=0 



(15) 



Note that the inclusion of the cos((n — A)7r/2) term is necessary in case n is an integer 
and n — A is even. In that case, sin((r;, — A)7r/2) = and the asymptotic behavior is 



described by the second term only. Eq. 15 is very useful as it gives a simple approximation 
to the Q-dependence of {j\{Q)) for large arguments. This is important as it allows to easily 
determine the maximal Q to be used in numerical simulations (or equivalently to estimate 
the systematic error introduced by considering only momentum transfers up to a certain 
maximum Q). 



For Q — i- the behavior of eq. 13 is more generally more complicatecP^, though. On the 
one hand, Q is nested deep inside several trigonometric functions, on the other hand many 
low-order terms cancel due to the unique form of the A'l^,B^^. In addition, evaluating eq. 



13 



numerically is also dangerous (because of the division by small numbers). Hence, for this 
case, evaluating eq. 12 for small Q directly is favorable. Provided that QRma.x <^ a/A + 3/2 
(where -Rmax is the supremum of radii for which the product iPi{R)iI)f{R) is non-negligible), 
this becomes 

poo 

NiNf I R''e-^^jx{QR)dR 
Jo 

NjNfQ^ 



(2A + 1)!! 
NiNFTin + X 



1)) Q' 



(16) 



C"+^+i (2A + 1)!! 

This shows that the usually used approximation eq. [9] is perfectly recovered for small 



Q. Furthermore, eq. 16 is important for the actual implementation in simulation software 
packages. As was noted before, inserting Q = into eq. [13] directly would produce a 



division-by-zero error. Hence, for Q ^ 0, eq. 16 should be used instead. 



III. HYDROGEN-LIKE ORBITALS 

Hydrogen-like orbitals are very similar to STO in so far as the same terms that appear 
for STOs appear in HLOs as well. There are four key differences, however: (a) in HLOs, 
the principal quantum number n is an integer, (b) the radial part of the wave function 
depends on the angular momentum quantum number /, (c) HLOs have an additional factor 
represented by (generalized) Laguerre polynomials, and it can thus be ensured that (d) 
HLOs with same /, but different n are orthogonal. 

In general, the radial part of the HLOs can be written as 

iVe-S(^)'e.,(^). (17) 

where the L'^^l_^{x) are the generahzed Laguerre polynomials, and = aorrie/ with the 
electron mass me, the reduced mass /i, and the Bohr radius acP'. Typically, ( is set to be 
equal to the (unscreened) nuclear charge Z. This does not influence the further calculation 
here, however, and so we use ( to indicate that screening may be included, and to preserve 
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the analogy to the STO results. The normalization constant is given by 



N 



2C V(^-/-l)! 



^na^J 2n{n + 1)\ 

With this model, the weighted radial wave function overlap eq. [8] becomes 
(jxiQ)) =NiNf / jx{QR)e 



' (I . Cf ' 



^2 f'KiB\' (KfRX 



di?, 



which can be rewritten as 



nj+np 

NjNf J2 R'e-''^''jxiQR)dR, 

b=lj+lp+2 

with C := { — — I — ^ ) and the coefficients Ph as defined in append: 



B 



(18) 



(19) 



(20) 



The integral has exactly the same form as for STOs, with the solution (see appendix |A 



{jxm 



NrN^ 



Q 



A+2 



nj+np 

E 

b=lj+lp+2 



Pb 



A+l 

2 



A\ sin ( {2k + 6 — A) arctan 



Q 



+ 



Bx,bVC + COS (2A; + 6 - A - 1) arctan 



(21) 



Since, contrary to the situation for STOs, b is an integer here, this can be simplified further. 
With 



sm 



cos 



n arctan 
n arctan 



^ [(C + 

(^2 + g2)n/2 

^ [(C + iQ)"] 

(^2 + Q2)n/2 



(22) 



one gets 



A+l 

nj+np 2 



PbQ 



2k 



QX+2 ^ ^ 1(2 



Q 



2\2k^b-X 



{A\m [(C + \Qf''^'-^\ + Bl,{e + Q^)3fJ [(C + iQ) 



2fc+fe-A-l 



]]• (23) 
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TABLE II. Weighted radial wave function overlap in the HLO model for monopole (A = 0), dipole 
(A = 1), and quadrupole (A = 2) transitions. For the definitions of the constants Nj , Np,n,(^ 
refer to the text. The first column shows the spectroscopic notation, the second shows the orbitals 
involved. 



K Is 

Is - 
L 2s 
2s - 
2s ■ 
2p 
2p 

2p 



2s A 
2p A 
3s A 
3p A 
3d A 
■ 3s A 
3p A 
A 

3d A 



NiNp 
NiNf 

(Q2+C2)2 

NiNfQ 

WW? 

NiNf 
(Q2+C2)2 

NiNfQ^ 

(Q2+C2)2 

NjNfQ 



2(Q^-30Cf 



2(Q^-3C^)(3C/+4Cf) 



12C + 

3a„(Q2+C2) 
128CC|. , 



64Cf(Q^-5C^)(Cf+3C/) 



16CCf(C'-Q')(2Cf+9Cj) 

+ 3a2(Q2+^2)2 



16C/Cf{Q^- 



3a3(Q2+^2)3 



+ 



9a2(Q2+^2)2 
64C7C|(Q2-7C2) 



64CCjC|(3Q^-5C^) 



24CC/ 



192CCjCf(C^ 



+ 



3a2(Q2+^2)2 



16C/Cf(Q'-5C^) _ 
32QC/Cf(Q^- 



32CC/Cf(3Q^-5C^) 
3aJ(Q2+C2)3 



+ 



128CC/CF 
a2(Q2+^2): 

64CCiC|^(5C"-3Q2) 



3a3(Q2+^2)3 



This is the second main result. It should be emphasized that eq. 23 has the form of a the 
quotient of two polynomials, 

(24) 



Q^+2((2 + Q2)n,+nF+l- 

As such, it is only marginally more complicated than the dipole approximation (which is 
simply a linear polynomial in Q), gives but an exact expression of the weighted radial wave 
function overlap in the framework of the HLO model. 

In table |TT| eq. 23 is evaluated for the most important transitions with A = 0, 1, 2, corre- 
sponding to monopole, dipole, and quadrupole transitions, respectively. 



As before, one can study the behavior of eq. 23 for the limit of Q oo and Q 0. For 



Q — 7- oo, eq. 21 reduces to 



ni+UF 



Pb 



b=lj+lF+2 ^ 



A+l 



sin ((&-A)|)X:(-l)' K. + i^y- 



A:=0 



c 

Q 



A+l 
2 



COS ((6 - A) I) [{2k + b - A)4, + (2fc + & - A - 1)5^ 



k=0 



■ (25) 



This exhibits exactly the same behavior as eq. 15 



For Q — 7- 0, it is best to start from eq. 20 Using eq. [oj (provided that QRmax ^ a/A + 3/2 
as before), one obtains 

N.N,^^ y (26) 
' '^(2A + 1)!! ^ C'+^+i ^ ^ 

in accordance with predictions. 



IV. DISCUSSION 



Eqs. 13 and p3^ are the main results of this work. Compared to the dipole approximation 



for small Q (eqs. [16] and 26), they offer a significant improvement. Unphysical cut-offs are 
no longer necessary to ensure they tend to zero properly for large Q. In addition, they 
are simple finite sums of polynomials and — in the case of STOs — elementary trigonometric 
functions. As such, they are easily implemented into existing simulation program^^^^ with 
very little effort and the increase in computation time is small. 

It must be emphasized that STOs are designed to have the same asymptotic behavior as 
HLOs, but no nodes. In the case of I = n — 1 (in which case the HLOs are also nodeless), 
they are identical. In all other cases, the missing nodes are problematic. This is particularly 
evident in the case of U = If which can be allowed in monopole or quadrupole transitions. In 
those cases, STOs give completely wrong results because the node-less radial wave functions 
are not orthogonal as they should be, as can be seen from fig. [TJ 

In fig. [2| examples of weighted radial wave function overlaps using both the STO and 



HLO models are shown in more detail. It is clearly evident that using the complete eqs. 13 



and 23 gives highly superior results than the approximations for small or large Q alone. 

The formulas even give comparable results to much more sophisticated calculations using 
a full crystalline environment of the atom. In fig. [3| a calculation using the HLO model is 
compared to WIEN2lP^ calculations using the TELNES.3 program. The excellent agreement 
can primarily be attributed to the fact that the initial state — which has a high probability 
density in close proximity to the nucleus — can be viewed as a filter on the final state. 
Consequently, the weighted radial wave function overlap is dominated by the shapes of the 
orbitals close to the nucleus, and crystal effects like bonding play only a secondary role. 




FIG. 1. Wavefunctions in the STO (top) and HLO (middle) models for 2p and 3p orbitals, as well 
as their overlap, together with the radial transition amplitude (bottom) for the quadrupole-allowed 
transition between those two states. Due to the missing node in the 3p wavefunction, the STO 

overlap is is shifted towards larger r and underestimated in strength. Consequently, (j2(Q)) is 
shifted towards smaller Q and severely underestimated in the STO model. 

V. CONCLUSION 

In this work, we have demonstrated that the transition probability in Coulomb scattering 
can be separated in an angular part and a radial part. The former is composed of well-known 
matrix elements of spherical harmonics. The latter can be expressed in simple independent- 
atom models such as the STO or the HLO model. Both yield simple algebraic expressions 
involving only finite sums of polynomials and (in the case of STOs) elementary trigonometric 
functions. Therefore, they are easy to implement into existing simulations programs without 
a large increase of the computation time, but with significant improvements in terms of 
accuracy, especially in the regime of medium momentum transfer. 

Moreover, some of the weighted radial wave function overlaps have one or more zeros. 
Hence, they are suppressed at the corresponding momentum transfers, even though they 
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FIG. 2. Weighted radial wave function overlaps {jx) for the example of an isolated Si atom. The 
top panel shows a 2p — )• 3d dipole transition calculated using the STO model (for Q > Q, the curve 
is magnified by a factor of 50) , the bottom panel shows a 2p — )• 3p quadrupole transition calculated 
using the HLO modelj^ The screening constants for both cases were taken from ref. [S The dashed 
lines show the asymptotic behavior. 

The different orders of magnitude are caused by the different screening constants for the different final 
states. For the 3d state, the nuclear charge is almost completely screened, resulting in a redistribution of 
the orbital's probability density to larger radii, and hence to a reduction of the overlap with the initial 
state, compared to the transition to a 3p state. For the 2p 3d transition, HLO and STO calculations 
are identical. 

are not forbidden by other selection rules. This can be exploited, e.g., for measuring faint 
signals from transitions that are normally hidden under a huge background from another 
transition with much higher transition probability. 

Finally, the formulas presented here can be exploited in the future to experimentally 
determine properties of wave functions in atoms, like the screening effects of other electrons. 
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Q (1/a.u.) 



FIG. 3. Comparison of the 2p — )• 3d dipole-transition calculated for crystalline Si using the HLO 
model (solid line) and WIEN2lil2l (dots). The screening constants for the HLO model were deter- 
mined by fitting to the WIEN2k data. For Q > 5, the data is magnified by a factor of 50. 
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Appendix A: Weighted integral of spherical Bessel functions 



In this section, we calculate 



(Al) 



For this, we use the (finite!) expansion of the spherical Bessel functions in terms of 
trigonometric function^ 



A+1 
2 



k=0 

sm(x) 



k\ 



'2 ^ 



^A-2fc+l \\-2k 

A+1 

2 



2k- 



COS X 



'2 ^ 



\-2k+2 \X - 2k + I 



k=0 



~u smlx) ~u cos(x 



^ rji.\-2k+l ~'~ >^ rji.\-2k+2 



where the binomial coefficients have to be understood in a generalized sense, i.e., 

T{n + 1) 



kj F(A; + l)F(n - A; + 1) 



(A2) 



(A3) 
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with the gamma function r(x). Since 



2 J {2k)\ ' 

the binomial coefficients can be simphfied to conventional factorials, yieldin 



B 



k 



{\-2k)\{2k)\ 
lf2'^k-\-i^(^2\-2k + 2)\ 



^ {\-k + l){\-2k + l)\{2k)V 
In particular, we have the special cases 

^°A = |^ = (2A-1)!! 5° = 



l)tA(A + l) 



> 2 



_ A+l „ A+1 



A+1 



A' =0 B,^ =(-1) 

and the general identities 

Xk+i ^ (A-2/c)(A-2fc-l) 

^ (A- A;)(2A-2A;- 1)(A; + 1)(2A; + 1) ^ 

(A-2fc + l)(A-2A:) 
= -^Tzr-, — —B^ 



k{2\ - 2k + - k){2k + I)' 



~, _ k{2X -2k + l) ~, 
- x-2k + l 

With this, we can writd^ 



E 

k=0 



R^e-'^^jxiQR)dR 

I {QRT^^'"^-^e-Q^'^'^hm{QR)d{QR) + 





A+1 



Q 







Qk 



Q 



n+l 



^ I {QR)''+^''~^-\-i^'^'^hos{QR)d{QR) 







A+1 

2 / /n2 



1 Q 



Al^^Q sin ^{2k + n - A) arctan ^ + 
BinVC + Q^ cos (^{2k + n-\-l) arctan (^^^ ^ 
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with 



T{2k + n- X)Al 



(A9) 



Appendix B: Polynomial coefficients of HLOs 



The polynomial part of the weighted radial wave function overlap using HLOs is given 



by 

\niaj ynpaj '''-''-^ \niaj ''^-'^-^ \nFaJ 
By virtue of the series expansion of the associated Laguerre polynomial^, 

''n + q\ 



j=0 



n- J J 3 



and the Cauchy product of (finite) series, this can be written as 



ni+np — li—lp — 2 k 
fc=0 j=0 



Up + Ip 

np — Ip — 1 ~ k + j 

ni+np 
b=li+lp+2 

E 

j=0 



+If+2, 



Tlj - Ij - 1 - j 



ripa^ 



jKk-jy. 



( 20 V' f 2Cp Y 



rij + li \ / np + lp 
ni-lj - I- J J \np + li + l- b + j 



which is of the form 



j\{b-k-lf-2-j)\ 



rij+np 



=ll+lp+2 



(Bl) 



(B2) 



(B3) 



(B4) 



lstefan.loeffler@tuwien.ac.atl 



^ H. Geiger and E. Marsden, |Proc. R. Soc. Lond. A 82, 495 (1909) 
^ E. Rutherford, Philosophical Magazine. Series 6 21, 669 (1911) 



14 



^ Of course, a scattering event can also produce entangled states. Since these can be decomposed 
into a linear combination of product states, the formalism presented here can be applied to 
those states as well. 

^ \i) and |/) can always be described in a plane wave basis. In the most general case, each is given 
by a coherent superposition of several plane wave components. Thus, the formulas derived in 
this article apply to each component of the general case as well. 

5 J. C. Slater, |Phys. Rev. 36, 57 (1930) 



^ E. Clementi and D. L. Raimondi, |The Journal of Chemical Physics 38, 2686 (1963) , 
''' R. F. Egerton, Electron Energy-Loss Spectroscopy in the Electron Microscope, 2nd ed. (Plenum, 
1996). 

^ M. Nelhiebel, Effects of crystal orientation and interferometry in electron energy loss spec- 
troscopy, Ph.D. thesis, cole Centrale Paris, Chatenay-Malabry (1999). 



^ M. Nelhiebel, P.-H. Louf, P. Schattschneider, P. Blaha, K. Schwarz, and B. Jouffrey, Phys. 



Rev. B 59, 12807 (1999) 



10 



P. Schattschneider, S. Rubino, C. Hebert, J. Rusz, J. Kunes, P. Novak, E. Carlino, M. Fabrizioli, 



G. Panaccione, and G. Rossi, [ Nat. 441, 4 86 (2006)) 

S. Loffler, I. Ennen, F. Tian, P. Schattschneider, and N. Jaouen, Ultramicroscopy 111, 1163 
1(2011) 



12 D. Essex, P. Nellist, and C. T. Whelan, ptramicroscopy 80, 183 (1999) 

It is possible to do the expansion on a case-by-case basis, though, such as for those cases listed 
in tab. |lj This yields the same results as derived from eq. [12] by using the approximation for 
small Q, of course. 

1^ For most practical applications, fj, ^ nie and hence = ao can be used. 
1^ S. Loffler and P. Schattschneider, jUltramicroscopy 110, 831 (2010) 



C. Koch, Determination of Core Structure Periodicity and Point Defect Density along Disloca- 
tions, Ph.D. thesis, Arizona State University (2002). 



1^ E. J. Kirkland, Advanced computing in electron microscopy (Plenum Press, 1998). 

1^ J. Verbeeck, P. Schattschneider, and A. Rosenauer, Ultramicroscopy 109, 350 (2009) 



19 



P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and J. Luitz, "Wien2k, an augmented plane 



wave + local orbitals program for calculating crystal properties," [(2001), (Karlheinz Schwarz, 



Techn. Universitat Wien, Austria). 

15 



M. Abramowitz and I. A. Stegun, eds., Handbook of mathematical functions (Dover, 1965). 
The integrals converge only if2A; + n — A>— 1 (for the sine integral) and 2k + n — X > 1 (for the 
cosine integral). Since n := rij+np, k > 0, and A < max(n/, np) = np without loss of generality, 
this translates to n/ > —1 for the sine integral, which is always fulfilled. Analogously, for the 
cosine integral we have k>l (since 5^ = 0), and hence the resulting inequality is identical. 



16 



